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Average entanglement of random pure states of an N x N composite system is analyzed. We compute the 
average value of the determinant D of the reduced state, which forms an entanglement monotone. Calculating 
higher moments of the determinant we characterize the probability distribution P(D). Similar results are ob- 
tained for the rescaled N th root of the determinant, called G-concurrence. We show that in the limit N — » oo 
this quantity becomes concentrated at a single point G* = 1/e. The position of the concentration point changes 
if one consider an arbitrary N x K bipartite system, in the joint limit N,K — » oo, K/N fixed. 
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O ■ I. INTRODUCTION 

fc: 

' Designing protocols of quantum information processing one usually deals with some particular initial states. One is then 
interested in describing the evolution of such a concrete quantum state and its properties in time. For instance, one studies the 
time dependence of the degree of quantum entanglement, which characterizes the non-classical correlations between subsystems 
and is treated as a crucial resource in the theory of quantum information [1]. 

As a reference point one may compare the degree of entanglement of the analyzed state with analogous properties of a typical, 
random state. Such random states are also of a direct physical interest since they arise under the action of a typical quantum 
^ chaotic system - see e.g. [2]. In this work we investigate mean values of certain measures of quantum entanglement, averaged 
' over the entire space of pure states of a Hilbert space of a given size. 

There exist several measures of quantum entanglement which do not increase under local operations and satisfy the required 
properties listed in [3, 4], but it is hardly possible to single out the "best" universal quantity. On the contrary, different en- 
tanglement measures occurred to be optimal for various tasks, so it is likely we will have to learn to live with quite a few of 
^ ■ them [5, 6]. 

Qh| The measures of quantum entanglement for a pure state of a bipartite system, \if/) e = 'Ha ® "Kb, rely on its Schmidt 
■ coefficients [7] equivalent to the spectrum A of the reduced system, p = Tib(\^){^\)- By construction the sum of all Schmidt 
coefficients equals unity, 2Zf=i A,- = 1, so just (N — 1) of them are independent. To quantify entanglement of a pure state one uses 
entanglement monotones [8], defined as quantities which do not increase under Local Operations and Classical Communication 
(the so-called LOCC operations). Entanglement of a pure state of a NxN system is therefore completely described by a suitable 
set of (N - 1) independent entanglement monotones. 

It is convenient to work with the ordered set of coefficients, Ai > A2 > ... > A^ > 0. The first example of such a set of 
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entanglement monotones found by Vidal consists of sums of k largest coefficients, Et := Yn=\ A; with k = 1, . . .,N — 1 [8]. 
Alternatively, one can use Renyi entropies of (N - 1) different orders. Another set of monotones may be constructed out of 
I symmetric polynomials of the Schmidt coefficients of order k = 2, . . . , N [9], 
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For large these polynomials become small, so it is of advantage to consider cognate quantities, r' k = (r^) 1 ^. Gour noted 
that taking the A^-th root of the polynomials does not spoil the monotonicity and proposed to used normalized quantities ri as 
alternative measures of quantum entanglement [10]. In particular he found unique properties of the last polynomial Tn, equal to 
the determinant of the reduced matrix, D = det p. Its rescaled AMh root, 

G:=ND UN , (1.1) 
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proportional to the geometric mean of all Schmidt coefficients, was called G-concurrence in [10], where its operational inter- 
pretation as a type of entanglement capacity was suggested. This quantity extended by the convex roof construction for mixed 
states, played a crucial role in demonstration of an asymmetry of quantum correlations [11] and was used to characterize the 
entanglement of assistance [12]. 

The aim of this work is to compute mean values and to describe probability distributions for the determinant D and its root 
G of random pure states of a bipartite system, generated with respect to the natural, unitary invariant measure on the space of 
pure states, also called Fubini-Study (FS) measure. Our analysis is performed for a bipartite system of an arbitrary size N, and 
in particular we treat in detail the interesting limiting case N — > oo. Although our study directly concerns bipartite systems, one 
may infer some statements valid also in the general case of multipartite systems. 

The paper is organized as follows. In section II we review a concept of a random pure state and describe certain probability 
measures in this set. Average values of the G-concurrence are computed in section III, while the subsequent section concerns 
with probability distribution of this measure of quantum entanglement. The paper is concluded with some final remarks while 
the discussion of the asymptotics of probability distributions is postponed to an appendix. 



II. RANDOM PURE STATES AND INDUCED MEASURES 



Consider a pure state of a bipartite N x K system represented in a product basis 
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The Schmidt coefficients A, coincide with the eigenvalues of a positive matrix pj^ = AA 1 , equal to the density matrix obtained 
by a partial trace on the TT-dimensional space. The matrix A needs not to be Hermitian, the only constraint is the trace condition, 
Tr AA' = 1. Furthermore, the natural unitarily invariant measure on the space of pure states corresponds to taking A as a matrix 
from the Ginibre ensemble [13]. Thus our problem consists in analyzing the distribution of determinants of random Wishart 
matrices AA 1 " normalized by fixing its trace. Schmidt coefficients's distributions are given by [14] 
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in which the cases of real or complex A are distinguished by the repulsion exponent /3 [15] being equal 1, respectively 2 and the 



normalization B^' K reads [13] 
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(2.1b) 



Formulae (2.1) describe a family of probability measures in the simplex of eigenvalues of a density matrix of size N. The integer 
number K, determining the size of the ancilla, can be treated as a free parameter. 

Another important probability measure in the space of mixed quantum states is induced by the Euclidean geometry and the 
Hilbert-Schmidt (HS) distance. Assuming that each ball of a certain radius contains the same volume, one arrives at the HS 
measure [16] 
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where the parameter B distinguishes as before between the real and the complex cases. The above normalization constant 
reads 



j- : = \ n 
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We observe that the distribution (2.2), normalization constants included, can be recasted into the form (2.1), provided that we 
choose K = N - 1 + 2//3, that is 



K 



\N for complex p^, (with /3 = 2) 
\N + 1 for real p N , (with/3 = 1) 



(2.3) 
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Using this observation, one can get a useful procedure for generating random density matrices distributed according to the 
HS-measure taking normalized Wishart matrices AA T , with A belonging to the Ginibre ensemble of Hermitian matrices of 
appropriate dimension. 

Aiming to derive the averaged moments needed in Section III, it is convenient to change variable in (2.1) by putting K = 
2a IP + N - 1 and obtaining: 



P^\A U ...,A N ):= C%*>6 
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In the above formula the real variable a can be used as a free parameter instead of the integer K. 



(2.4b) 



III. AVERAGE MOMENTS OF G-CONCURRENCE 



In this Section we are going to compute averages over an ensemble of random density matrices distributed according to the 
HS-measure, which is induced by the Euclidean geometry. This corresponds to fixing the size K of the ancilla according to (2.3), 
depending on whether the real or the complex case is concerned. 

Denoting the eigenvalues of the density matrix p N by {A;}, the moments of the determinants D(Ai , . . . , A N ) = JT^i ^-j rea d 
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The product of Heaviside step functions, present in the definition (2.2a) of Pj®, allows us to extend the domain of integration on 
the entire axis. The integrand of (3.1) coincides with the factor present in the right hand side of equation (2.4a), provided that 
the parameter a is set there to 1 + M. Using this the integral (3.1) can be computed from (2.4b), and reads 
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For sake of clarity, from now on the sub- and super-script (J3 = 2) and (/} = 1) will be often replaced by C, respectively R. 
Making use of equation (1.1), one obtains the moments of the G-concurrence by imposing a — 1 +M/N in the ratios /C^\ 
rescaled by a factor N M . Thus we get now 
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(3.3) 



in FIG. 1 the mean values (G^) and variance cr 2 N = (G^) - (G^) are represented as a function of N for both complex and 
real cases. 



IV. PROBABILITY DISTRIBUTION (G) 

This Section is devoted to the study of probability distributions. We shall start with the simplest problem of determining the 
distribution of the determinant D of a 2 x 2 density matrix p2 distributed according to the (HS)-measure. In this case an explicit 
solution is easily obtained by integrating the Dirac delta 6(D - A1A2) over the distribution pj^(Ai, A2) of (2.4), that is 

dA, I dA 2 5(A, +A 2 - l)6>(Ai)0(A 2 ) |A, - A 2 f 6(D - A, A 2 ) • 

CO %J —OO 
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FIG. 1 : Average of G-concurrence for (a) complex and (b) real random mixed states of a N x (N + 2 - 0) system distributed according to 
the HS measure. The average is computed by means of equation (3.3); error bars represent the variance of P^(G). Dashed line represent the 
asymptote G* = 1/e, whose explanation is given in Section V. 



It is a very simple distribution since (A] - A 2 ) 2 = (A] + A 2 ) 2 - AD = 1 - 4D. Thus 

D e 



\p^{D) = 6 Vl -4D 



Pf(D) = 4 



(4.1) 



The G-concurrence distribution pf\G) can be computed either by integrating 6(G - 2 VA1A2) over P^g(Ai, A 2 ), or simply 
using the latter result (4.1) together with pf{G) AG = pf{D) AD; in both cases (see FIG. 2) 

{f G) = 3G ^ , Ge[0( l]. (4.2) 
\Pf(G) =2G 

Note that, due to Ai + A 2 = 1, (only) for the case N = 2 the G-concurrence given by (1.1) reduces to the standard concur- 
rence [17], C = -^2(l - Trp^). Thus formula (4.2) for the complex case coincides with the distribution of concurrence P(C) 
obtained in [13]. 

For higher N we will construct the distribution P^iD) from all moments (D^)n given by equation (3.2); indeed 



dDD M P%\D) = J o dxe 



with D = e 4 and AD = -Ax e x , and so we can obtain P^{D) by inverse Laplace transform or inverse Mellin transform as 
integral along the imaginary M-axis: 

^ )(D)= irs D_(i+M)<jD ^ (4 - 3) 

ns (4.3) and (3.2) allow us fc 

as well by using 



Although equations (4.3) and (3.2) allow us to compute the P^\D) probabilities, the cognate quantities P^\G) can be determined 



Pf(G)AG = Pf(D)AD; 



by taking from (1.1) the explicit expression for AD /AG, one can indeed get the simple expression 

G • Pf(G) = N ■ D • Pf{D) ■ (4.4) 

From now on formulae and figures will be given indifferently for both G and D distribution, being clear their mutual relation. In 
particular the D-distribution is more indicated in showing details of calculation, for its simpler form, whereas the G-distribution 
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FIG. 2: G-concurrence's distributions Pfi (G) are compared for different N in the case of (a) complex and (b) real random pure states. The 
distributions are obtained by performing numerically [18] the inverse Laplace transform of equation (4.3). Dashed vertical line centered in 
G+ = 1/e denotes the position of the Dirac-delta corresponding to Pf^iG), as it is shown in section V. 



better shows features in the pictures, for its domain being independent of N. Important is the asymptotic behavior of the Gamma 
function for large argument (Stirling's formula) 



for |z| — > oo and arg(z) < n. This implies the asymptotic behavior of (3.2) for large \M\: 
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As a consequence the integral (4.3) converges and moreover it vanishes if x < NlogN or D > (l/N) N , because in that case we 
can close the contour in (4.3) in the right M-halfplane according to the Jordan's Lemma [19]. Physically this means that there 
are no density matrices with determinants greater than the one with maximal entropy. 

In the rest of this section we will give the asymptotic behavior of distributions P ( ^(D) for the two edges of the domain, that 
is D — > and D — > (1/N) N . The details of calculation, together with the explicit N-dependence of all coefficients listed here in 
the following, are collected in Appendix A. 

In particular, when very close to the completely mixed state, that is D (1/N) N , we have the result (see FIG. 3) 
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Moreover, using (4.4) together with 



we simply find 
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For the other part of the spectrum, that is for very small D, the probability P^(D) can be expanded in a power series with some 
logarithmic corrections, as follows: 

P%(D) - Z% + X% ■ D log D + X% ■ D + V% ■ D 2 (log D) 2 + V% ■ D 2 (log D) + V% ■ D 2 + O (D 3 (log D) 3 ) (4.7) 

In particular, coefficients Z^j, X®, X% are computed in appendix A for all Af > 3, whereas for V®, and we limit ourself to 
explicitly solve the case N = 3 (the case N = 2 is simply given by formula (4. 1)). 

The situation is similar when we do consider, in the same region of the domain, the probability P^(D), corresponding to small 
determinants of reduced N x N real density matrices HS-distributed. The expansion is still a power series (plus logarithmic 
corrections) but the exponents are now semi-integer, according to the mechanism described in Appendix A, thus the probability 
reads: 

P%(D) + T* • D l 2 + X^ ■ DlogD + X^ ■ D + ■ logD + ■ + O (z) 2 (log D) 2 ) 



V. CONCENTRATION OF G-CONCURRENCE FOR LARGE SYSTEM SIZE 



Iterating the recursion relation for the Gamma function T(n + 1) = n T(n), we can recast expression (3.3) of the M- moment 
of the G-concurrence of complex random pure state as 
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FIG. 4: In panel (a), formula (4.2) is compared with a 100 bins histogram of 10 6 G-concurrence of 2 x 2 complex density matrices distributed 
according to the HS measure. The other panels (b) (c) and (d) shows histograms (for different AO together with the distribution of G- 
concurrence obtained by inverse Laplace transforming as in equation (4.3) (plotted in solid lines). The left asymptote given by eq. (4.7), 

computed up to O (D), is also plotted in dashed line for comparison; in panel (b) we also add the contribution given by V®, Vfj, coefficients, 
using a dotted line. 



so that finally 

(Gj?) »e- M • (5.1) 



For the analogue moments of G-concurrence of real random pure state, some technicality requires that the sequence of odd and 
even N has to be analyzed separately, although it is not hard to prove that the limit is the same. For that reason, we will simply 
illustrate the case N = 2p, p e IN, for which (3.3) gives 
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Putting all factors together we arrive at the general result (compare with (5.1)) 

G(M):= iTJ G A = e " • (5 - 2) 

The above expression, valid for both (3 e {1, 2), is useful to derive the limiting distribution 

P W (G) := lim Pf{G) ■ 

We see from (5.2) that its average is 1/e = 0. 367 879 441 . . . and its variance is 0; such behavior can be recognized in FIG 1. 
Moreover, by fixing G = e~*, one can see that G(M) of (5.2) is nothing but the Laplace transform of the function 

j](x) := e^P^e - *) 



so that, by inverse Laplace transforming, we obtain 

e -*pW( e -*) = G pW((j) = <J(-log(G) - 1) 
Rewriting the argument of the Dirac delta we finally arrive at 



P W (G) = S(G - e" 1 ) 



(5.3) 



In other words, we have shown that for large systems the G-concurrence of random states is localized arbitrarily close to the 
averaged value. 

A similar concentration effect has recently been quantified [20] for bipartite N x K systems. In particular the Von Neumann 
entropy of the reduced density matrix of the first subsystem concentrates around the entropy of the maximally mixed state, 
S (1 IN) = log N, if we let the dimension K of the auxiliary subsystem to go to infinity faster than N. When K — N, so that 
the induced distribution coincides with the Hilbert-Schmidt distribution, and N — > oo, then von Neumann entropy concentrate 
around log N - 1 /2 [20, 21]. Remarkably, G-concurrence displays a similar concentration effect; moreover, we are in position 
to prove the convergence of its distribution to a Dirac delta centered at a non trivial value 1/e. 

The determinants and G-concurrence may be also averaged in the general case of asymmetric induced measure (2. 1). Consider 
an interesting case K > N. As for the HS-distribution discussed in Section II the expectation value and the higher moments may 
be expressed as a ratio of normalization constants (2.1b) and (2.4b). For instance, the moments read 
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Let us now study a particular case of the induced measure, for which we consider bipartite systems of arbitrarily large dimension, 
with the only constraint that the ratio between the size K of the ancilla and the size N of the principal subsystem are fixed and 
greater than one. Let this ratio be expressed by the rational number q = till\ , with the C\ and (2 integers; this means that we are 
considering systems with N - J(\ < K - Jli. 

With the same tools used in computing (5.2), one can let J go to infinity and obtain 

G(M):=hm(G-)^ 2 =X-- , V/>e{1.2} , (5.5) 

with 

X q := I(^p . q >l ■ (5.6) 
The limiting distribution P 9 (G), can be earned as before and reads 

Pf (G) := lim Pf (G) = 6(G - X q ) , 

for the complex as well as for the real case. Although the accumulation point X q is not defined for the case q = 1 (that is the 
case in which states in the principal system are HS-distributed), we find however lim 9 ^i X q = 1/e, confirming our previous 
result (5.3). Moreover such values represent an infimum for X q , whereas it attains the supremum on the other part of the domain, 
that is for q — > 00. Such case correspond an extremely large environment, for which Xca = 1, that is in turn the G-concurrence of 
the completely mixed state. Thus we find another evidence that large environment concentrates reduced density matrices around 
the maximally mixed states [20]. 

VI. CONCLUDING REMARKS 

The generalized G-concurrence is likely to be the first measure of pure state entanglement for which one could find not 
only the mean value over the set of random pure states, but also compute explicitly all moments and describe its probability 
distribution, deriving an analytic expression in the large N limit. This offers for our work various potential applications. On one 
hand, analyzing a concrete quantum state and its entanglement we may check, to what extend its properties are non typical. In 
practice this can be done by a comparison of its G-concurrence G with the mean value (G), and by comparing its deviation from 
the average, |G - (G)|, with the root of the variance of the distribution. 

On the other hand, if one needs a quantum state of some particular properties, one may estimate how difficult it is to obtain 
such a state at random. For instance, looking for a state of a large degree of entanglement, with concurrence greater than a 
given value G, one can make use of the derived probability distribution by integrating it from G to unity in order to evaluate the 
probability to generate the desired state by a fully uncontrolled, chaotic quantum evolution. 

Although in this work we have concentrated our attention on pure states of bipartite systems, the averages obtained for the 
asymmetric induced measures (2. 1) with K > N may be easily applied for the more general, multipartite case. Consider a system 
containing n qudits (particles described in a <f-dimensional Hilbert space). This system may be divided by an arbitrary bi-partite 
splitting into m and (n - m) particles, and one can study entanglement between both subsystems - see e.g. [22]. The partial 
trace over m qudits is equivalent to the partial trace performed over a single ancilla of size K = d m , so setting size of the system 
N = d"~ m one may read out the average concurrence from eq. (5.4). In particular, if n is even and we put m = n/2 + k, then 
the ratio q = K/N is equal to d 2k and in the asymptotic limit n — > 00 the concurrence concentrates around the mean (5.6) which 
depends only on the asymmetry k of the splitting. 

Our research may also be considered as a contribution to the random matrix theory: we have found the distribution of the 
determinants of random Wishart matrices AA 1 ", normalized by fixing their trace. Furthermore, the analysis of the distribution 
of G-concurrence in the limit of large system sizes provides an illustrative example of the geometric concentration effect, since 
in high dimensions the distribution of the determinant is well localized around the mean value. This observation can also be 
related to the central limit theorem applied to logarithms of the eigenvalues of a density matrix, the sum of which is equal to the 
logarithm of the determinant. 

It is a pleasure to thank P. Hayden and P. Horodecki for stimulating discussions. This work was financed by the 
SFB/Transregio-12 project financed by DFG. We also acknowledge support provided by the EU research project COCOS 
and the grant 1 P03B 042 26 of Polish Ministry of Science and Information Technology. 
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APPENDIX A: COEFFICIENTS OF ASYMPTOTIC EXPANSIONS OF PROBABILITY 
1. Right asymptote of P^(D): proof of equation (4.6) 

The starting point is integral (4.3). Since all the poles of the integrand are in the left half-plane (see it in (3.2)), the contour 
integration along the imaginary axis can be modified into the one along the right asymptotic half-plane, that is on a very 
large semicircle connecting -zoo to +;oo; this allow us to use the Stirling's formula for replacing {D^)n with (M,N) (see 
formula (4.5)) in the integrand of (4.3). Of course we made an approximation, but we know that the formula we ended up 
matches the correct result (P^'(D) = for D > ( 1 /N) N ) in the point ( 1 /N) N , so that such approximation would hold close to that 
point. Now we observe that (M, N) has poles only in M = 0, so that our contour of integration can be modified provided 
that we do not cross the origin, and we do so obtaining 

P V\ D) = C^L D-d+M) D S (M N) = i£. f dM e M(-lo S D-Nlo g N) M -^l 

N J 7 2m m 2mD J y 

where y is now the contour that, starting from -zoo get close to the negative real axis on the asymptotic left-lower quarter-plane, 
winds around R~ U {0} in the counterclockwise direction, and then approaches +z'oo on the asymptotic left-upper quarter-plane. 
But now we apply once more Jordan's Lemma and we remove the asymptotic semi-circle from y. After rescaling M — > -M/s, 
with the latter defined by s — -logD - NlogN and close to 1, we arrive at the well known Hankel's contour integral for the 
inverse of the Gamma function (I/O [23], that leads to (4.6) and gives the asymptotic behavior for D — » (1/N) N . 



2. Left asymptote of P®(D) for complex random pure states 



Now let us consider the behavior of P^(D) at the lower edge of the spectrum D — > 0. In that case one can close the 
integral (4.3) in the left halfplane obtaining contributions from all the poles of the Gamma functions in (D^)jv ( see (3-2)). Such 
poles are located at each of the negative integers M = -1,-2, -3, . . .; fortunately there is the factor £)-( 1+M ) such that we obtain 
a series in powers of D. Because of the multiple Gamma functions in (3.2), most of the poles are degenerate and the general 
feature (for an arbitrary large N) is that the pole in —( is of order t: due to this fact the D-powers in the expansion get in general 
a logarithmic correction. The first pole at Mi = - 1 is non degenerate and yields 



T(N 2 ) 



T(N 2 - N) F(N) N 

Including the next order-2 pole {Mi = -2) contribution we find the asymptotic expansion for D 

P%(D) - Z% + X% ■ DlogD + 5£ • D 



with 
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N T(N 2 - 2N) T(N) T(N - 1) 
X% = X%(N + Ni//(N 2 -IN)- A- 2i/r(l) - (N - 2)iff(N - 2)) 

Here ifr(x) is the Digamma function [24], or polygamma function of order 0, with 



(Al) 



^(1) = -J, <A(«) = -y + V -, for n > 1 . 

U k 



(A2) 



Note that the Euler constant y cancels everywhere. By adding the next order-3 pole (Mj = -3) contribution one gets in 

general the terms in (4.7) corresponding to the V®, and coefficients, although the latter are in general rather complicated, 
involving polygamma function of order higher than 0. This is not the case when N = 3, for which a cancelation makes M3 = -3 
a pole of order 2, and the coefficients read : 



=0 



6-7! = 30240 



yV = 9-7! =45 360 
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3. Left asymptote of P^(D) for real random pure states 

We will apply the same reasoning of the previous case, just now differing for the fact that, when /3 = 1, the & pole Mi of 
the integrand of (4.3) is in general, for arbitrarily large N, its corresponding order is given by [^y-J, where [xj means the 
larger integer not exceeding x. In particular, the firsts two poles M\--\ and M 2 = -\ are non degenerate and yield [25] 



7 R _ 



2 N-l T{ N^N ) 
r(^) T(N) 



and 



jN-l 



r(^) 



r(^v) r(*H) r(Ar-i) 



Including the next two 2-order poles contributions (M 3 = -2 and M 4 = -|) we determine, for N > 3 



yR 



r |«^wj r(A , )r(A ,_ 2) 

:^{jv + ^(^)-8-§^)-2^(1)-^^)-^^)) 



(A3) 



and for N > 4 



W R 



W R 



3 ^/fl^N j r(^i) r(Af-i) r(Jv-3) 
<(Ar + ^(^)-f-f^(i)-2^ ( i)-^^(^)-^^(^)) 



(A4) 



where we made use once more of the i^-digamma function [26] of (A2). The case N = 3 constitutes an exception for X R and 
Wf 's coefficients, because of the lowering of the order of M3 and M4 poles; moreover, for the latter pole, the same happens also 
for N - A. All these coefficients need separate calculations and read 



Xf 



12-5! 



= 4-5! 



W 



■K 



2 s • 8! 



and 



= W: 
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